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HEAT TRANSFER CHARACTERISTICS OF AN EMERGENT STRAND 

By William Emile Simon, Larry Claude Witte,* 
and Pat Garner Hedgcoxe* 

Lyndon B. Johnson Space Center 

SUMMARY 


The results of an analysis of the heat transfer characteristics of a hot strand 
emerging into a surrounding cooler fluid are presented . The analysis is for a stable 
strand of constant efflux velocity, with a constant heat transfer coefficient on the lead- 
ing surface and strand sides. A dimensionless governing equation is derived and solved 
by Laplace transform methods. The equation provides a description of the variation of 
temperature within the strand with axial distance and time. Following a discussion of 
computational aspects of the solution, generalized results are presented for a wide 
range of parameters. A discussion of the relationship of the generalized results and 
experimental observations is given. 


INTRODUCTION 


When a hot emergent strand of an incompressible fluid or semisolid material pen- 
etrates a cool liquid bath, heat is lost rapidly from the sides and leading edge of the 
strand. If the' strand is a hot molten material, explosive vapor formation may occur, 
caused by fragmentation of the molten material (ref. 1). Conversely, the strand may 
be of a metal or a fiber in its plastic deformation region, being extruded and quenched 
simultaneously. 

In the area of nuclear reactor safety, an over -temperature excursion within a 
reactor can cause overheating of the fuel elements, possibly to the point of becoming 
molten, while also generating fission gases and increasing the pressure level inside 
the fuel elements. When this happens, the temperature of the cladding that surrounds 
the fuel element rises, ultimately causing a cladding rupture that can result in the re- 
lease of a high-temperature molten jet stream into the surrounding coolant. Film boil- 
ing begins immediately on the exposed material and continues until the stream cools 
sufficiently for the system to enter the transition and subsequent nucleate boiling re- 
gimes. The violent nature of transition and nucleate boiling quickly overcomes the 
surface tension forces holding the molten metal together and tears the molten mass 
into smaller globules. The process sets off a thermal chain reaction because, as more 
exposed surface area is provided as a result of disintegration, the heat transfer 
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increases, causing further disintegration by nucleate boiling. The fragmentation proc- 
ess continues until, in some cases, the molten fluid is reduced to a very fine powder. 

If the time to transition boiling can be predicted, information can be obtained on the 
amount of energy released up to this point and on the delay times involved in the vapor 
explosion process. 

The transportation of liquefied natural gas (LNG) by large oceangoing vessels is 
another application of the processes analyzed in this report. A leak in the hull of an 
LNG vessel could allow seawater to penetrate the LNG. Because of the large tempera- 
ture difference between the two fluids, the LNG would begin to boil off the sides of the 
seawater stream, resulting in a phenomenon analogous to the analysis of this report. 

Another application is in the metal-forming industry, in which it is often neces- 
sary to form metal bars by means of an extrusion or, perhaps, a rolling process. A 
closely related situation concerns a metal wire being drawn from a die, or a monofila- 
ment fiber being drawn from a spinneret (melt spinning). As indicated in reference 2, 
only in the recent past has melt spinning been the subject of fundamental research. If 
the cooling behavior of a wire (or billet) being extruded could be predicted, it might be 
possible to combine the metal-forming and the heat-treating processes into one. 

The applications just described are the ones that provided the motivation for this 
study. In these applications, a fundamental understanding of basic heat transfer within 
the strand is essential. Many parameters govern the thermal performance of the strand. 
In this report, a mathematical model is developed and the governing equation is derived 
to describe the heat transfer characteristics of the strand. The governing equation is 
then solved for the temperature within the strand as a function of axial distance and 
time . Finally, a computer program is used to vary the parameters governing the heat 
transfer process and to determine the degree of influence of each. (The applicability 
of this analysis to emerging liquid strands or streams is very limited because, by as- 
suming a stable strand geometry, the complex fluid dynamic interactions between the 
strand and the surrounding cooler fluid are neglected.) 

For this report, a stable, one -dimensional, circular strand of constant efflux 
velocity is assumed, with a constant (average) heat transfer coefficient on the sides and 
leading surface of the strand. However, the heat transfer coefficient on the end of the 
strand need not be equal to that of the sides. For convenience, the cylindrical body is 
referred to as a strand, whether it is a wire, billet, fiber, molten jet, or other form. 

The approach to the heat transfer solution is similar to that of the classical fin 
heat transfer analysis (ref. 3), but the problem is complicated by the fact that the strand 
is in motion, while at the same time it is always physically connected to the fixed reser- 
voir, resulting in a constantly changing geometry. Therefore, a departure from classi- 
cal techniques is needed to arrive at a solution to this particular problem for the short 
times of interest in some applications (e.g. , the vapor explosion problem). For large 
times, the solution is analogous to the well-known semi-infinite rod problem. 

As an aid to the reader, where necessary the original units of measure have been 
converted to the equivalent value in the Systeme International d’Unites (SI). The SI units 
are written first, and the original units are written parenthetically thereafter. 


2 



SYMBOLS 


Biot number, hr/k 

Biot number at the end, h g r/k 

specific heat of strand material 

/ 2 

Fourier modulus, ar/r 
time to steady state 

heat transfer coefficient on strand sides 
thermal conductivity of strand material 
Peclet number, rU/a 
Prandtl number, c^p/k 

heat generated within elemental volume 

heat entering elemental volume 

heat lost from strand sides 

heat leaving elemental volume 

heat stored in elemental volume 

local heat transfer rate per unit area 
strand radius 
Reynolds number, rUp/p 
transform variable 
quadratic roots 

temperature 
reservoir temperature 


coolant temperature 



U strand efflux velocity 

u Heaviside unit function 

x coordinate position, distance from the entrance 

z function, defined in equation (B9) 

a thermal diffusivity of strand material 

y dimensional parameter, 2h/kr 

rj dimensionless distance, ij/r 

r)* steady-state distance 

0 dimensionless temperature 

M absolute viscosity of strand material 

1 position coordinate, distance from the leading edge 

p density of strand material 

t time 


FORMULATION OF THE ANALYTI CAL MODEL 


To begin the analysis, consideration first is given to the flow of a round axisym- 
metric strand in the region of stable film boiling with constant efflux of the molten ma- 
terial or strand. As opposed to the analysis of Schlichting (ref. 4) for a jet that mixes 
with the surrounding fluid, in this analysis a stable jet or strand is assumed, with no 
mixing of the strand material and the coolant. Neglecting leading-edge deformation 
(ref. 5 is suggested as a starting point for analyzing the effects of velocity retardation 
and deformation), the physical process can then be modeled by a blunt-edged strand of 
radius r penetrating the coolant as shown in figure 1. A vapor film forms at the lead- 
ing surface and continues toward the wall, completely engulfing the strand sides. At 
any given position, the vapor film thickness will grow with time, always being thinnest 
at the leading edge. Because the film boiling process greatly reduces skin friction at 
the strand surface (refs. 6 and 7), velocity gradients within the molten material are 
neglected. In addition, because of the high thermal conductivity of liquid metals and the 
small radii of practical interest, radial temperature gradients within the strand are ne- 
glected, In summary, the following assumptions are made concerning the problem: 
a stable axisymmetric strand at constant velocity, radial temperature gradients ne- 
glected, and constant heat transfer coefficient. 

When a hot or molten strand emerges into a region of cool liquid, three factors 
are of primary importance: the fluid-dynamic characteristics of the strand system; 
the boiling of the coolant; and the cooling and solidification of the molten material. For 
the case of a fiber or strand being drawn through a coolant medium, boiling may or may 
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Wall 

Regions; A ■ molten material 

B s vapor film (coolant vapor} 

C s coolant liquid 

Figure 1.- Emergent strand penetrating 
surrounding coolant. 
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not be involved, and the cooling process is 
simpler than that of the molten material be- 
cause no solidification is involved. 

For the molten strand problem, it is 
assumed that, when the strand first emerges 
from the opening, stable film boiling with 
strong radiation effects occurs (ref. 8). As 
the strand travels downstream, cooling as 
it goes, it eventually reaches a temperature 
at which stable film boiling exists without 
appreciable radiation effects. Further cool- 
ing results in a region of partial nucleate 
boiling and unstable film boiling. Next 
comes the nucleate boiling regime; if the 
strand does not disintegrate, pure convec- 
tion will follow. The various heat transfer 
regimes are shown in figure 2 (ref. 9). 
Experimental work concerning vapor explo- 
sion indicates that under certain^conditions 
the molten material will disintegrate vio- 
lently somewhere in the region of nucleate 
boiling. The disintegration and subsequent 
solidification process is very complicated. 
Disintegration produces a rapid increase 
in heat transfer area that increases the heat 
flux, thus causing a self-propagating effect 
and hastening the heat transfer process even 
further. In addition to the various heat 
transfer regimes, the velocity regions of 
laminar and turbulent flow are also impor- 
tant for both the molten material or fiber 
and the cooler liquid (liquid and vapor 
regions). 


Figure 2.- Boiling of water at 373° K 
(212° F) on an electrically heated 
platinum wire, showing various 
boiling regimes. 


To determine the heat transfer coef- 
ficient governing the removal of heat from 
the surface of the strand, the equations of 
heat, mass, and momentum transfer must 
be accompanied by appropriate boundary 
conditions for each heat transfer and veloc- 


ity region, and further boundary conditions connecting the various regions must be sat- 
isfied. If the heat transfer coefficient is known beforehand or if an average coefficient 
is used, the analysis is concentrated on the heat transfer within the hot material, which 
is the case for this report. More information on the calculation of the heat transfer 
coefficient is given in appendix A. (Since the heat transfer coefficient is treated para- 
metrically in the main body of this report, this appendix is therefore an independent 
entity and can be considered entirely separate from the report. The purpose of appen- 
dix A is to merely suggest one possible approach to the extremely difficult problem of 
computing the heat transfer coefficient.) 
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To derive the governing differential equation of the strand system, consider first 
the elemental cross section of the strand shown in figure 3. In the general case, h is 
the overall coefficient of heat transfer for convection, radiation, and film boiling. The 
quantities q i and are those of heat transfer into and out of the elemental volume. 

Heat stored within the element is denoted by q g ; q £ is the heat lost off the sides be- 
cause of convection and so forth. It is assumed that the temperature of the surrounding 
coolant T m remains constant. ° 


It is convenient to translate the coor- 
dinate axes to a moving coordinate system 
at the leading edge of the strand. This 
translation is accomplished by means of the 
— t*. — transformation 



{ ■ Ur - x 


(1) 


x + dx 


where £ is the position coordinate, dis- 
tance from the leading edge; U is jet ef- 
flux velocity; r is time; and x is the 
coordinate position, distance from the en- 
trance. At all times, £ > 0. The moving 
coordinate system is also shown in figure 1. 
An elemental heat balance including heat 
conducted into and out of the element, heat lost from the sides of the element, and net 
heat accumulation by the element yields the basic equation 


Figure 3.- Elemental volume of strand 
showing energy balance. 


4 

3? 2 


- re 


= ±de 

a 3 t 


( 2 ) 


where y = 2h/kr, a dimensional parameter; a is the thermal diffusivity of the strand 
material; and the dimensionless temperature q is defined as 


T - T w 

o(t) s - - (3) 

0 °o 

where T is the temperature of the strand material and T is the reservoir 
temperature. ° 
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For equation (2) to apply in both regions (outside the wall on the strand surface 
and inside the reservoir), the Heaviside unit function u (ref. 10) is introduced. 



( 4 ) 


If the actual temperature gradient at the entrance is modified so that q = e Q = 
at £ = Ur , then multiplying the heat transfer coefficient y in equation (2) by 

u[-^(r - |f)] results in the following modification of equation (2): 
r 

In this formulation, when £ < Ur , heat is transferred from the sides at the rate y; 
when £ > Ur , no heat is allowed to flow from the sides of the strand within the reser- 
voir (dotted lines in fig. 1). 

In seeking a closed-form analytical solution, it is apparent that in the transformed 
coordinate system, in which the position coordinate £ is always greater than or equal 
to zero, the governing equation is particularly suited to the Laplace transform method 
of solution. However, solution of equation (5) by the Laplace transform method leads 
to a more complicated equation, and another approach is therefore taken. 

The energy balance equation rewritten in its most general form is q. - q Q + q gen = 

q s + q £ ( where the heat g eneration term ^en iS Zer ° outside the wal1 )* However, 

inside the wall, because the term q £ dictates that heat still is being lost from the sides 

of the strand, an equal and opposite amount of heat must be generated within the strand 
to replenish the heat lost. This make-up quantity of heat can be expressed as 

q =2?rrhd£^T -T 
H gen \ o 
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because the strand is essentially at T 0 
for this physical model then becomes 


inside the reservoir. The governing equation 


( 7 ) 

where g * (T q - Tj/(T 0 - Tj = 1. 

In reality, equation (7) is more physically correct than equation (5) because it allows 
for a temperature gradient at the entrance. Outside the entrance, u(r - ~) = 1, and 
equation (7) is identical to equation (5). However, within the entrance, where 
u[-s (t - ^y)] = 0 , the two equations differ by the term y(g - g ); and, as £ becomes 

r v ° 

large, T must approach T , in which case convection from the side vanishes, as 
expected. 

Equation (7) may be put in dimensionless form by the introduction of the following 

variables: 77 = £/r, Fo = ar / r 2 , Bi = hr/k, and Pe = rU /a = (rUp/ju)(c p ju/k) = RePr 

where 77 is the dimensionless distance, Fo is the Fourier modulus, Bi is the Biot 
number, k is the thermal conductivity of the strand material, Pe is the Peclet num- 
ber, p is the density of the strand material, ju is the viscosity of the strand material, 
c p is the specific heat of the strand material, Re is the Reynolds number, and Pr is 

the Prandtl number. Then, equation (7) becomes 

0- 2Ble =^S- 2Bi [ 1 - u ( Fo -^)] (8) 



If equation ( 8 ) is applied to a solid wire or billet, the Peclet number representation is 
meaningless. However, Pe is basically a nondimensionalized velocity rU/o'. 

Two boundary conditions and an initial condition are imposed as follows. Bound- 
ary conditions: Fo > 0, 7 ? = 0; || = Bi e e , where Bi e is the Biot number at the end; 

other boundary condition is 7 } - e is bounded. Initial condition: Fo = 0, 77 = 0; 

0 = 0 Q = 1. With this formulation of the analytical model, attention is now directed 

toward obtaining a solution for this system that will accurately predict the physical 
phenomenon. 
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METHOD OF SOLUTION 


If the Laplace transform of both sides of equation (8) is taken with respect to Fo 
(denoting the transform £{e| by 9 = 0(77, s), where s is the transform variable), the 
result is a linear second-order ordinary differential equation in 0 and r ) . 

- (2B1 + e)e = -1 - ?fl[l - exp(- §£)] 
dr? 

The complementary solution can be written as 

e c = A(s)exp(?7>j2Bi + s) + B(s)exp (-7?>j2Bi + s) 
where A(s) and B(s) are to be determined from boundary conditions. The particular 

solution is 


( 9 ) 


( 10 ) 


0p = C^s) + C 2 (s)exp(-srj/Pe) 


( 11 ) * 


The evaluation of the constants from the boundary conditions yields the general solution 
of the transformed equation 


e (r?,s) = exp^T^J^BT 


- 2B1 (l 

, Bi ) 
1 e' 

3 e + s ) 



- 

/ s > 

_W 

2 

■ 

• (2Bi + s) 

(Bi g + >/2Bi + s) 


Bi e 


1 x - 

2Ei exp(- 


{ + c 

It 

r 2 ■ 

s(Bi e +>j2Bi + sj 

1 ^ 1 


(A) - (2Bits) J 


( 12 ) 


This solution was verified by taking the appropriate derivatives and by substituting them 
and equation ( 12 ) into equation (8). The inversion of equation ( 12 ) provides the general 
solution for the strand system. 
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The inversion details are omitted for the sake of brevity. The details, although 
representing a considerable portion of the effort in solving this problem, are generally 
straightforward and are given in reference 11. The solution is lengthy and involves ex- 
ponential and complementary error functions. It is written as 


e(7?,Fo) = 1 - u^Fo - - 2BiPe exp(-2BiFo)^A 1 + B^PeD^ ^ exp^Fo) 

UBip 7 ? + Bi p Fojerfc/ — - — + Bi \[Fo\ 

V 6 6 J \2 yjFZ eX ) 


L( Bi e + \Z i T) 


exp (- 


Bi 2 - s 
. e 1 


2BiPe exp(-2BiFo)^B 1 + Bi g PeEjJ 


) 5 ex P( s 2 Fo ) 

1 

(Bi + Jsj) 


\ e > 2J 


• exp (- , 'fi) er,c (^ - 'J 5 ? 5 ) * “ p ( , 'Pi) 


Bi JF^ 


+ 2BiPe 2 E^ exp 


2BiPe 2 D 1 exp - 2Bij ^Fo - 

it - 2Bi > (fo - &)] 


(-’SP7) 


where 


A l = 


Pe(pe 2 + 8Bi) 


(14) 


B l = - A l 


(15) 
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D„ = 


1 

r - 

^(pe 2 + 8Bi) 2 

Pe 2 + Pe(pe 2 + 8Bi) _ 


(16) 


E l = 


^ (pe 2 + 8Bi) Pe 2 - Pe (pe 2 + 8Bi) 


Tt 

2 


(17) 


and where and s 2 are roots of a quadratic given by 


12 = ^4Bi +Pe 2 ) ± [(4Bi +Pe 2 ) 2 - 16Bi 2 ] 


(18) 


The ultimate check of correctness for the general solution was performed success- 
fully by taking the appropriate partial derivatives from equation (13) and substituting 
them into the governing differential equation (8). 


APPLICATION OF THE SOLUTION 


The solution, equation (13), yields e as a function of rj and Fo, with Bi, Bi g , 
and Pe as parameters. The Biot numbers Bi and Bi e represent the ratio of internal 

to external resistance to heat transfer and must be kept quite low if the neglect of radial 
temperature gradients is to hold. For this investigation, Biot numbers were limited to 
0. 1 or less. 

Equation (13) is applicable for all times; however, as time (Fo) becomes large, 
arguments in the solution become correspondingly large. If these terms are separated 
from the overall solution and are designated for convenience e 2 , then 


0 2 


= 2BiPeD 1 exp 


• erfc 





+ Pe exp 




(19) 
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The remainder of equation (13) would then be called e j, so that e = e ^ + ©g- This 

grouping of terms eliminates complications caused by machine limitations at large 
times. Generally, as Fo - 00 , the exponential functions involving Fo become larger. 

It can be shown that for long times, approaching steady state, e 2 becomes of negli- 
gible significance. A brief demonstration of this is given in appendix B. 

Figure 4 shows the end temperature ratio of the strand e(0,Fo) for a Biot num- 
-3 

ber of 10 and for various Peclet numbers. The lower the Peclet number, the longer 
the time required for the end to cool to essentially the coolant temperature. Initially, 
this fact seems contrary to intuition; however, the effect of axial conduction is being 
demonstrated here. Low Peclet numbers correspond to strands that move away from 
the reservoir slowly. Axial conduction is effective in these cases, and the strand cools 
more slowly. Figure 4 also shows that high Peclet numbers (Pe > 1) have little effect. 
This situation corresponds to the case of a more rapidly moving strand where axial 
conduction does not contribute so heavily. 

For small metal jets or strands, the Peclet number is fairly high except for very 
slow moving strands. Figures 5 and 6 show the cooling results for the leading edge for 
two Biot numbers smaller than that of figure 4. When the leading edge approaches the 
temperature of the coolant (e = 0), no further cooling of the strand past that point occurs. 
Essentially, figures 4, 5, and 6 yield the Fo at the time when the strand approaches a 
steady -state condition. At this point, the temperature profiles in the strand undergo 
no further change. 



strand: Bi = 10 - *\ 



Figure 5. - Cooling of the end of the 
strand: Bi = 10 
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0(0, FO) 



Fo 


Figure 6.- Cooling of the end of the 
strand: Bi = 10 \ 



Pe 


Figure 7.- Time required to reach 
steady state for low- Biot- number 
strands . 


The time to steady state Fo* is 
determined by 

0 ( 11 , Fo, Bi, Pe) = e(0, Fo*, Bi, Pe) = 0 

( 20 ) 

In actuality, e only approaches zero as- 
ymptotically; therefore, Fo* was deter- 
mined as the time required for e to reach 
9 = 0.01. Figure 7 shows the steady-state 
time plotted with Bi as a parameter for a 
range of Pe that will encompass most situ- 
ations of practical interest. The figure 
shows that Fo* decreases with increasing 
Bi but increases with decreasing Pe for 
Pe < 10. For Pe > 10, the Fo* is essen- 
tially constant, indicating the same effect as 
shown in figure 4 but in a slightly different 
manner; that is, relatively large Pe has 
a negligibly small effect on the cooling rate 
of the strand. Clearly, the Bi could be a 
function of Pe, but that interdependence was 
not included in this investigation. 

However, the Pe does have a bearing 
on the length of the strand when it reaches 
a steady -state condition. The product 
PeFo* gives the steady-state distance rj* 
as the strand tip approaches the coolant 
temperature. Figure 8 shows the variation 
in steady- state length for a wide range of 
Pe. This figure illustrates that, if Bi and 
Pe are known for a strand, the nondimen- 
sional 77 * can be found easily. 

Figures 9, 10, and 11 show how the . 
temperature profiles develop as time pro- 

2 

ceeds for strands with Pe = 10 . The dis- 
tance scale is normalized versus the length 
of the strand at a particular Fo, For ex- 
ample, in figure 9, the plots for Fo of 

5 x 10 1 , 5 x 10 2 , and 10 3 are for Fo < Fo*; 

4 

and the curve for Fo = 10 is for Fo > Fo*. 
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Figure 8.- Length of strand required 
for steady state. 



Figure 10.- Development of tempera- Figure 11.- Development of tempera- 
ture profile: Bi = 10" 2 , Pe = 10 2 . ture profile: Bi = 10“\ Pe = 10 2 . 


The test of any theoretical result is a comparison with experimental results. 
Aciemo et al. (ref. 2) have obtained temperature profiles for polyethylene fibers at 
t0 « 73 o K / 180 ° to 200 ° C ) bein s pulled from an extruder into stagnant air at 
293 to 298 K (20° to 25° C). The conditions for this test do not exactly match those 
of the theory developed in this report; for example, the strand was being extended and 
e diameter decreased as the extension took place. Therefore, the velocity and diam- 
eter of the strand are functions of the distance from the spinneret. Another difference 
is possibly due to the variation in thermal conductivity in the radial direction as a result 



of radial temperature gradients. Even with these differences, axial temperature distri- 
butions computed from the theory of this report show an excellent comparison in shape 
with those measured by Acierno, indicating that the theory is correct, at least in its 
dependence on basic parameters of the system. 

Bradley (ref. 12) undertook a companion experimental study to this theoretical 
study. He studied the cooling characteristics of a 0. 159-cm-diameter (1/16-inch- 
diameter) molten metal jet in distilled water. The effects of drag caused the jet to be 
slowed as it penetrated into the liquid, resulting in considerable waviness in the jet. 
Direct temperature measurements along the jet could not be obtained, but in some 
cases behavior was observed that compares well with the theoretical results of this 
analysis. Bradley noted that solidification of the jet occurred at a fixed axial location 
from the injection nozzle, after the jet had penetrated the cooling bath considerably. 
Under similar conditions, the theory of this investigation predicts that the temperature 
at which solidification would begin is reached at a constant axial distance from the 
reservoir. 

/ 

Data that could be. compared directly to the theory were not found. Complete veri- 
fication of the theoretical formulation must await such results. Meanwhile, work is 
being conducted to obtain data by means of an apparatus involving the passage of a solid 
wire through a cooling region. 


CONCLUDING REMARKS 


: A mathematical model was developed to describe the heat transfer characteristics 

of a hot jet or strand emerging into a surrounding coolant. A closed-form analytical 
solution of the governing partial differential equation was obtained for a circular strand 
of constant efflux velocity and constant (average) heat transfer coefficient on its sides. 
Using this model, the various parameters governing the thermal behavior of the strand 
were investigated. Although a limited comparison with experimental data was made, 
more experimental data are needed with which to compare the theory. 

The complexity of determining an adequate heat transfer coefficient on the strand 
was illustrated by a simplified analysis consisting of an extension of previous work to 
include several effects that would be present in most applications of interest. In the 
model, it was necessary to use an average heat transfer coefficient. The effects of 
variations in this coefficient were investigated parametrically. For a specific physical 
application, an initial average heat transfer coefficient can be computed, and this aver- 
age coefficient can then be used in the model to determine the thermal performance of 
the strand. With this information, the computation for the heat transfer coefficient can 
be refined. This procedure can be repeated as many times as necessary until a suffi- 
ciently accurate average overall heat transfer coefficient for the strand system is 
obtained. 

The principal limitations of the analysis are summarized as follows. For the 
case of a freely moving strand, the assumption of a stable axisymmetric strand at con- 
stant velocity does not provide for velocity retardation or deformation of strand mate- 
rial as it penetrates the coolant. This limitation is not applicable in the case of a fiber 
or strand being pulled through the coolant at a constant velocity . These effects perhaps 
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would not be significant for the small times of interest for fast-cooling small strands 
because stable film boiling would soon give way to transition and subsequent nucleate 
boiling, which is beyond the realm of this report. However, if the temperature at 
which transition occurs is known, then the solution will provide a prediction for the 
time required for a strand to experience transition. The neglection of a radial tempera- 
ture gradient within the strand has the effect of limiting the application to materials of 
high thermal conductivity or to strands of small radius. The assumption of a constant 
average heat transfer coefficient necessarily implies that an accurate value of this co- 
efficient be computed, and this computation is the one that most severely limits the 
analysis for any given application. 

Therefore, areas of further study should concentrate on striving to eliminate the 
previously mentioned limitations so that the scope of application of this effort may be 
widened. The inclusion of a radial temperature gradient within the strand leads to the 
consideration of Bessel functions. It would allow for more accurate analysis of strands 
with large radii. Further work is needed in the area of determining the heat transfer 
coefficient for this geometry. The uniform- temperature assumption presents a signif- 
icant limitation, and much work is needed here, both experimental and analytical. In 
particular, hydrodynamic aspects and their effect on the heat transfer process should 
be investigated in greater detail. The effects of velocity retardation and deformation 
should be studied. Perhaps such a study could lead to deformation predictions and to 
the determination of a time-varying heat transfer coefficient at the leading edge. 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, December 10, 1973 
981-10-00-00-72 
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APPENDIX A 


v 

DETERMINATION OF HEAT TRANSFER COEFFICIENT 


In this appendix, a simplified analytical approach directed at calculating an aver- 
age heat transfer coefficient for forced-convection film boiling on the strand sides is 
presented. Although some of the same symbols that were used in the text will be used 
in this appendix, many of the meanings are different. 


SYMBOLS 


Cp specific heat at constant pressure 

F dimensionless liquid stream function, defined in equation (A14) 

f dimensionless vapor stream function, defined in equation (All) 

h overall local heat transfer coefficient, q”/(T w - T^) 

hjg latent heat of vaporization 

h local radiation heat transfer coefficient, defined in equation (A22) 

r 

k thermal conductivity of vapor 

m mass flow rate at interface 

Nu overall local Nusselt number, hx/k 

Nu radiation Nusselt number 

r 

Pr Prandtl number, c u/k 

r 

a" conduction heat transfer 

c 

a" radiation heat transfer 

^ r 

Re Reynolds number, Ux/V L 

T static temperature 
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T sat coolant saturation temperature 

/ 

wall temperature 

free-stream temperature 

U free-stream velocity 

u velocity component in x direction 

v velocity component in y direction 

x position coordinate 

y position coordinate 

£ parameter for subcooled boiling, defined in equation (A43) 

6 vapor film thickness 

e L thermal emissivity of liquid coolant 

€ w thermal emissivity of strand material 

V vapor similarity variable, defined in equation (A10) 
fjg dimensionless vapor film thickness 

e dimensionless temperature for vapor, defined in equation (A12) 

© L dimensionless temperature for liquid, defined in equation (A33) 

k thermal diffusivity 

X* modified latent heat, defined in equation (A23) 

P absolute viscosity 

v kinematic viscosity 

k liquid similarity variable, defined in equation (A13) 

P density 

1 P stream function for vapor 

stream function for liquid 
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HEAT TRANSFER COEFFICIENT ■ 

| 

Althou'gh a constant (average) heat transfer coefficient (which is treated paramet- 
rically) is used in obtaining the solution of the governing heat transfer equation (18) of 
the strand system, the overall coefficient of heat transfer on the sides of the jet is, in 
general, a complex function of the thermodynamics and hydrodynamics of the flow sys- 
tem with its associated boundary- layer interactions. Other factors that complicate the 
analysis are radiation effects and subcooled boiling. 

An analysis to determine heat transfer and friction characteristics during forced- 
convection film boiling on a flat plate, neglecting radiation and subcooling effects, was 
performed by Cess and Sparrow (ref. 7), and their approach is used in the analysis of 
this appendix; however, unlike the former study, this analysis includes the effects of 
radiation and subcooled boiling. In applying the flat-plate results to the sides of the jet, 
the effects of curvature are ignored. The work of Okabe (ref, 13) shows that the incom- 
pressible laminar boundary- layer flow along a two-dimensional flat plate is substantially 
similar to that of axial flow along a circular cylinder near the leading edge, where the 
ratio of boundary-layer thickness to strand radius is small. Conceptually, the cylindri- 
cal surface of the jet is approximated by incremental surface elements such as sec- 
tion AA in figure A- 1(a). Figure A- 1(b) illustrates the growth of the liquid and vapor 
boundary layers beginning at the leading edge of the surface element, and also shows 
the two-dimensional rectangular coordinate system used in the boundary- layer analysis. 


So that the governing equations will be analytically tractable, the following as- 


sumptions are needed: 



(a) Oblique view of emergent strand, 
illustrating applicability of flat- 
plate analysis to an axial sur- 
face element. 


1. A stable vapor film exists, with 
all generated vapor remaining in the film. 
Also, no interfacial waviness or film in- 
stabilities are present. 

2. The flow is steady, with negligible 
gravitational forces. 



(b) Analytical model for boundary- layer 
analysis, showing free-stream 
velocity U relative to emerging 
flat-plate element of strand. 

transfer coefficient on strand sides. 


Figure A-l.- Physical model for calculating heat 
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3. The plate surface temperature is uniform. Therefore, the determination of 
the heat transfer coefficient must be based on an average strand surface temperature, 
which is not accurately known. This implies an iterative process for obtaining a refined 
value of the overall heat transfer coefficient. 

4. A constant-property vapor is assumed, or a variable-property vapor charac- 
terized by pm = constant, pk = constant, and Pr = constant (where p is density, 4 
is absolute viscosity, k is the thermal conductivity of vapor, and Pr is the Prandtl 
number). 

5. Inertia forces and energy convection within the vapor film are neglected, and 
a linear velocity profile is assumed in the vapor layer. 

6. The density -viscosity product in the vapor layer is very much smaller than 
that in the liquid layer; that is, (pp) « (pm)^. * 

7. Regarding radiation effects, the vapor layer is assumed to be a nonpartici- 
pating gas; that is, absorption and emission within the vapor layer are neglected. This 
assumption is based on the work of Sparrow (ref. 14) for film boiling in free convection 
on a vertical flat plate, in which it is concluded that the effects of a radiatively partici- 
pating vapor on heat transfer are negligible within the parameter range investigated, 
which includes steam at pressures at least as great as 10 atmospheres. 

Because the derivation of the equations is quite lengthy, only a general discussion of 
the analysis is presented here. More detailed information on the method used is avail- 
able in the cited references. 

The analysis is begun by first considering the coolant fluid to be a saturated liquid. 
The effects of subcooling are incorporated later. In laminar flow, the governing equa- 
tions are those of continuity, momentum, and energy. They are expressed as 


9u 9v 
ax + ay 


0 


(Al) 


0u 3u 0 u 

U -r— + V -r— — V X 

dx 9y - 2 

0y 


0T 3T 0 2 T 

U -r— + V -t — = K if 

3x 3y 0 y 2 


(A2) 


(A3) 


In this appendix, variables without subscripts refer to the vapor, and the sub- 
script L is used for the liquid. 
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where u is the velocity of fluid in the x direction, v is the velocity of fluid in the 
y direction, x and y are position coordinates, v is the kinematic viscosity, T is 
the static temperature, and k is the thermal diffusivity. Although equations (Al), 

(A2), and (A3) apply to both the liquid and vapor, equation (A3) is not used for the liquid 
coolant in the preliminary formulation because the temperature in the liquid is assumed 
to be everywhere constant and equal to the saturation temperature. As written here, 
the governing equations represent a constant-property vapor. However, as shown in 
reference 7, for a variable -property vapor adhering to the relations PM = constant, 
pk = constant, and Pr = constant, the final heat transfer and skin friction results are 
identical to those predicted by the constant-property analysis. 

For continuity at the liquid-vapor interface (y = 5, where 5 is the vapor film 
thickness), the tangential velocities, tangential shears, and mass flows that cross the 
interface are equated for the liquid and vapor. 


Tangential velocity: u^ = u 


(A4) 


Tangential shear: ^(ly) = ^(^y) 


(A5) 


Mass flow crossing interface: P L |u ^ = ^ “ v ) (A6) 


For temperature continuity at the interface, the vapor temperature at the interface must 
be equal to the liquid temperature; that is, 


T = T at y = 5 

OO J 


(A7) 


where T^ is the free-stream temperature. 

Other boundary conditions needed are the conventional ones that apply at the wall and 
in the free stream. 


At y = o, u = v = o and T = T 


w 


(A8) 


As y — 00 , u — U 


(A9) 


where T is the wall temperature, and U is the free-stream velocity, 
w 
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Then, as in reference 7, boundary -layer similarity transformations are intro- 
duced, first for the vapor layer in which 


111 

DO hi 

it 5 ! 

(A10) 

f iv) - 

(All) 

y JT'UX 


T - T oo 


0 (^]) rp - T 

(A12) 


w 


where rj is the vapor similarity variable, f is the dimensionless vapor stream func- 
tion, is the stream function, and e is the dimensionless temperature. Then, for 
the liquid layer, 




y- 5 P0~ 

2 ^ L x 


(A13) 


F(£)s L . ... (A14) 

where £ is the liquid similarity variable, F is the dimensionless liquid stream func- 
tion, and \f/ 1 is the stream function. 

The similarity parameters are used to transform the governing partial differential 
equations into the following ordinary differential equations (ref. 15). 


f»"+ff" = 0 (A15) 

e" + (Pr)f 9 * = 0 (A16) 

F m + FF" = 0 (A17) 


where the primes denote differentiation with respect to tj in equations (A15) and (A16), 
and with respect to £ in equation (A17). 
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Next, transformation of the boundary conditions, equations (A4) to (A9), results 
in the following conditions, first at the heated surface 


f (0) = f’{0) = 0 and e(0)=l (A18) 


then at the liquid-vapor interface 


F(0) 


F”(0) 




(A19) 


where t}^ is the dimensionless vapor film thickness, and finally in the free stream 


F* — 2 as £ — °° 


(A20) 


Equations (A15) to (A20) form the complete set of governing equations for the trans- 
formed problem. Knowing the fluid properties, the three unknowns f, F, and e can 
be obtained by solution of equations (A15) to (A17) if something is known about This 

information is obtained by writing an energy balance at the interface. At this point, a 
departure from the approach of Cess and Sparrow (ref. 7) is necessary in order to in- 
clude radiation effects. When the sum of the conduction and radiation terms at the 
interface is equated to the heat required to evaporate the liquid, the energy balance 
becomes 


, 3T ./, 
-k -r— + h [ 
3y r \ 


T - T . 
w sat 


)■ 


mX* 


(A21) 


where h^ is the local radiation heat transfer coefficient, T ^ is the coolant satura- 
tion temperature, m is the mass flow rate at the interface, and X* is the modified 
latent heat. In equation (A21), 



(A22) 
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and 


X* B 



0.84 
1 + 


c (T - T V 
pV w saty 

h, • Pr 


(A23) 


where is the thermal emissivity of the liquid coolant, .e is the thermal emissivity 
of the strand material, h^ is the latent heat of evaporation, and is the specific 
heat at constant pressure. 


The radiation heat transfer coefficient (eq. (A22)) takes the stated form with the 
assumption that the vapor layer is very thin; and because of this, the radiant interchange 
between a surface element and the liquid occurs between locations that are essentially op- 
posite one another. For this reason, the nonparallelism of the surface and the inter- 
face will not affect the local radiant heat transfer, and, if curvature of the strand sides 
is neglected, then the configuration factor for two parallel planes can be used. 


The A* as defined in equation (A23) is used in equation (A21) because it has been 
shown (ref. 14) that the effects of the vapor inertia forces and of superheating essen- 
tially can be eliminated from the problem by defining such a modified latent heat of 
vaporization in terms of the actual latent heat h^, the temperature difference 


T - T , and the vapor properties, 
w sax 


fg’ 


Transforming the right side of equation (A6) for m into the new variables and 
evaluating the gradient at the interface, equation (A21) becomes 



(A24) 


where Re is the Reynolds number, and Nu is the overall local Nusselt number. The 
parameter ?jg is thus replaced by the dimensionless parameter c^AT/A*, which is a 

function only of physical properties and of the temperature difference AT = T - T ,. 

w sat 

According to Cess and Sparrow (ref. 7), certain simplifications can be made at 
this point that permit an easy yet accurate solution for the range of parameters en- 
countered in the applications of interest for this investigation. The first simplification 
is to neglect inertia forces and energy convection within the vapor film. This results 
in simplification of equations (A15) and (A16) to 


f ,M = 0 and e ” = 0 


(A25) 
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the solutions of which, together with the boundary conditions at the interface and plate 
surface, allow the determination of f(r? ) and e '(f? Q ), from which equation (A24) 
becomes 0 ° 


Equation (A26) reduces to Cess and Sparrow's (ref. 7) equation (18) when radiation ef- 
fects are negligible, with the exception that the modified latent heat A* is used instead 

of h. . 
fg 

The velocity problem in the liquid is identical to that solved by Cess and Sparrow 
(ref. 7) with the aid of assumption 6 mentioned previously; therefore, it need not be re- 
peated here. The method of solution relies on the adaptation of the Emmons-Leigh 
table of solutions (ref. 16) for the problem of blowing or suction normal to a flat plate. 
(See reference 7.) Several values for the derivatives F'(0) and F"(0) are shown in 
table A- 1 for use in equation (A26) and subsequently. 


"c at" 
p 

W) L ' 

_A*(Pr)_ 



2 


2Nu 


‘ PM 

• W) 


f ' 2 <°) me1 l/2 
fTW R 


f 


f^(0) (Re) l/2 


F T (0) 


(A26) 


TABLE A-I. * VALUES OF DERIVATIVES AND RELATED FUNCTIONS 
USED IN THE VELOCITY PROBLEM FOR THE LIQUID 


F '( D ) 

F M (0) 

F ,2 (0) 

F "(0) 

F "(0) 

fW 

0.000 

1.328 

0 


.491 

1.218 

. 198 

2.48 

. m 

1. 135 

.393 

1.70 

.787 

1.067 

.580 

1.36 

.877 

1.011 

.761 

1. 15 

1.005 

.926 

1.091 

.921 

1.117 

. 842 

3.482 

.754 

1.175 

.799 

1.728 

.680 

1.20 

.787. 

1.830 

.656 

1.30 

.700 

2.414 

.539 

1.40 

.612 

3.203 

.437 

1.50 

.519 

4.335 

.346 

1.60 

.423 

6.052 

. 264 

1.70 

.322 

8.975 

. 189 

1.80 

.218 

14. 863 

. 121 

1.90 

.111 

32.523 

.0584 

1.92 

. 0904 

40. 779 

.0471 

1.94 

. 0678 ♦ I 

55. 510 

. 0349 

1.96 

.0451 

85,180 

.0230 . 

1.98 

.0226 

173.469 

.0114 

2.00 ' 

.00 

00 

.00 
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The local heat transfer from the strand sides is that caused by conduction (q" c ) 
through the vapor film and by radiation (q” ) from the surface, and the local heat trans- 
fer rate per unit area (q M ) can be expressed as 


q" = q" 


+ q 



(pm) l 1 1 F „( 0 ) 

. pm J fW 


+ h r AT 


(A27) 


where the conduction term is that derived by Cess and Sparrow, and the radiation coef- 
ficient is as defined in equation (A22). The overall local heat transfer coefficient h is 
defined as 



(A28) 


Therefore, division of both sides of equation (A27) by AT yields the desired expres- 
sion for h. 


h = 



+ h 


r 


(A29) 


Using the definitions from the nomenclature for the combined radiation-conduction 
Nusselt number, the radiation Nusselt number Nu r , and the Reynolds number, an ex- 
pression is obtained for the Nusselt number caused by combined radiation-conduction. 


Nu 


(Re) 


172 



_ 1 F"(0) . Nu r / M \ 


(A30) 


Then, if equation (A26) is rearranged into the same dimensionless parameters as equa- 
tion (A30) 


c AT 
_P 

"<pm) l 

r*(Pr) 

PM . 


1 F ,2 (0) 

2 F»(Q) 


'2(Nu) r 

(M > 

F"(0)‘ 

W 751 

kJ 

1 + TW 


(A31) 


26 



Equations (A30) and (A31) can then be used in conjunction with table A-I to produce the , 
plot of figure A- 2, in which the influence of radiation is indicated. When radiation ef- 
fects are negligible, the bottom curve shown is identical to that obtained by Cess and 
Sparrow. 



Figure A-2.- Effect of radiation on 
overall heat transfer. 


The analysis will now be extended to 
include subcooled boiling. Where previously 
it had been assumed that the coolant liquid 
was saturated throughout (that is, = T gat ), 

it must now be specified that < T gat , 

but at y = 5 for any x, T(x, 6) = T gat . 

An additional equation, the energy equation 
for the liquid, is needed. 


u 


9T 

9x 


+ v 


3T 

3y 



(A32) 


By defining the dimensionless parameter 


0 




(A33) 


where e T is the dimensionless temperature, and by recalling the previous parameters 
L 

for the liquid layer 



(A13) 


F(0 


*_L 


(A14) 


on taking appropriate derivatives, equation (A32) is transformed into the following ordi- 
nary differential equation 


d^e T de T 

£ + (Pr) L F(e)- 5 | 1 = 0 (|>0) (A34) 

d4 
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or 


e L M + ( p r) L Fe L r = 0 (£>0) 


(A35) 


where the primes denote differentiation with respect to £. The necessary boundary 
conditions are for any x 


e L (0) = 1 


(A36) 


and 


lime (y) = 0 (A37) 

y-*oo 


Cess and Sparrow (ref. 17) have developed special relations for the temperature 
gradient at the edge of the boundary layer for three liquid Prandtl number regions. 


e L f (0) 

* L '«» 


e L^°) = 


= -21 


(Pr) T 


-F ,f (0) 

2 - F’(0) 


(low liquid Prandtl number) 


(liquid Prandtl number unity) 




2(Pr) L F’(0) 

7T 


(high liquid Prandtl number) 


(A38) 

(A39) 


(A40) 


If an energy balance is written at the liquid-vapor interface (y = 6 or £ = 0), it can be 
reasoned that the energy convected by the vapor and radiated from the jet surface goes 
into coolant evaporation and energy convected by the liquid at the interface; thus 

- k (f ) + h r( T w - T sat) * - k L (f) L < A «> 
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If the gradients in the liquid and vapor layers are evaluated at the interface in terms of 
the transform variables, equation (A41) becomes 


c AT 
P 


F ,2 (0) 


1 

2 F" 


|_A*(Pr)J ~ 2^Nuj /fl ~ 
(Re) 1 ' 2 v L/ 


m 


F"(0) 

fW 


+ /3 0 r ’(0) 


(A42) 


where /3 is a parameter for subcooled boiling and is defined as 





(A43) 


where 


AT s 


sat 


- T. 


(A44) 


Equation (A42) is used in conjunction with equations (A38), (A39) or (A40), and equa- 
tion (A30). Equation (A30) is still valid because heat is transferred from the strand 
surface both by vapor convection and by radiation, irrespective of subcooled boiling. 
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APPENDIX B 


SOLUTION SIMPLIFICATION FOR LARGE TIMES 


It can be easily shown that 




Then equation (19) of the text can be written, using equation (Bl), as 

« 2 = 2BiPeD 1 exp[( 5l - 2 b1)(fo-^)] [u(Fo-i)]pe 


- + Bi Pe 
\D, e , 


1 ' Bi e + >FI, 


1 \1 


2VFo 


-> 1 * 7 ° 


It also can be shown that 


2BiPe iuf + B1 e Pe 

i-J L - 2BiPe 


and because 


erfc(-x) = 2 - erfc(x) 


equation (B2) may be written as 


9 2 ’ 2BiPe 2 Dj exp[(Sj - 2Bi)( F o - £)] jj erfc^i 


2VFo, 


U ( F “-^)- 1 ] 
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In the region of interest, that is, where r\ < PeFo, 


u(Fo - 77 /Pe) -1 = 0 


(B6) 


and presents no problem. The remaining part of equation (B5) as Fo — 00 now must 
be examined. First, the square of the complementary error function argument is con- 
sidered; that is, 


/^sTFo - — \ 2 - 

V 1 2 W 


where equation (Bl) is used again. 



It is possible to write 



(B7) 


(B8) 


If 


z 


I {{ s 7° 


- 2-1 

2VF0/ 


(B9) 


then equation (B5) may be written as 


9 2 = 2BiPe D 1 exp 


2BiFo + 




(BIO) 


As Fo — *>, z — 00 . Reference 18 shows that 



if 


< exp 


2 

z + 2/ 


(z^)erfc(z) s -^L 
v 17 



(Bll) 
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2 

Clearly, as z — 00 , the bounding expressions for exp{z )erfc(z) approach zero. Con- 
sequently, for long times (large Fo), e 2 - 0 and e = e 1 only. For cases in which 

the jet or strand (in this case, probably a bar, wire, or fiber) is issuing continuously 
into the cooling region, only this portion of the solution need be used. 
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